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Theoretical study of systematics of neutron scattering cross sections on various materials for 
neutron energies up to several hundred MeV are of practical importance. In this paper, we analysed 
various cross sections of neutron-nucleus (n-N) scattering for several systems in the energy range 
of 50-250 MeV, predicted by the optical model using Koning-Delaroche potentials. We propose an 
empirical approach to successfully predict the energy dependence of total, shape elastic, reaction 
cross sections and zero degree scattering angular distributions. We demonstrate that owing to two 
conditions, only two out of these four cross sections need to be fitted empirically. Further, we 
modified the Wick's approximation for the zero degree angular distributions and use this approach 
for estimating various cross sections for any nucleus from Aluminum to Lead. 
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In recent times, the concept of an accelerator driven sub- 
critical (ADS) system is drawing worldwide attention 
PI Q- In this ADS system, neutrons are produced by 
bombarding a heavy element target with a high energy 
proton beam of typically above l.OGeV with a current 
of > 10mA PJ. Such a system serves a dual purpose 
of energy multiplication and waste incineration. In this 
context, theoretical study of systematics of neutron scat- 
tering cross sections on various materials for neutron en- 
ergies up to several hundred MeV are of practical impor- 
tance. In this paper, we propose an empirical approach to 
successfully predict the energy dependence of total cross 
sections (cr tot ), shape elastic cross sections (cr se ), reaction 
cross sections (a reac ) and zero degree scattering angular 
distributions (oo or a(9 = 0)), of the neutron-nucleus (n- 
N) scattering. The total cross sections are usually fitted 
by using the Ramsauer model. The nuclear Ramsauer 
model was first proposed by Lawson in the year 1953 as 
a simple means to understand the energy dependence of 
total cross sections of neutron nucleus scattering. In or- 
der to appreciate this model, it is necessary to discuss 
briefly the optical model (OM) description of neutron 
scattering. In the OM approach, complex optical model 
potentials (OMP) are used and the Schrodinger's equa- 
tion is solved to obtain the scattering amplitude. The 
real part of the OMP describes the scattering and the 
imaginary part results in attenuation or absorption of the 
incident wave. This absorption gives an estimate of the 
optical model reaction cross section. The calculations are 
usually performed using partial wave expansion method 
and the phase shifts (rj£ = a£e l ^ e ) are determined. These 
complex phase sifts are strongly angular momentum and 
energy dependent for a given set of potentials. In terms 
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of the phase shifts and the wave number (A = h/y2mE), 
various cross sections are given below. 
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Extensive study of the optical model fits of scattering 
cross sections on various nuclei over wide energy range 
have been made by several groups. This is owing to the 
excellent data base of neutron cross sections available in 
the energy range up to 600 MeV @, 0L@ . The most re- 
cent work by Koning and Delarosche 6] presents a very 
exhaustive search for OMP parameters that fit the data 
very well up to 200 MeV. In the remaining part of our 
work, we use the OM code SCAT2 with Koning and 
Delaroche potentials. We treat this as the "experimental 
data" for shape elastic, reaction, total scattering cross 
sections with spherical targets. We made a phenomeno- 
logical Ramsauer model analysis of this data to derive 
systematics of Ramsauer model parameters and utilize it 
for further predictions. As mentioned before, the nuclear 
Ramsauer model || provides a simple means to param- 
eterise the energy dependence of neutron nucleus total 
scattering cross sections. This model assumes that the 
scattering phase shifts are independent of angular mo- 
mentum (I) as given in Eq.(5) (77 = ae 1 ^) , in contrast 
to the predictions of the optical model given in Eq.(l). 
Further, it was proposed that the ^-independent phase 
shift varies slowly with energy. The model thus implies 
that the nucleus is seen as a right circular cylinder by 
the neutrons incident along the axis of symmetry. Ow- 
ing to this unphysical picture, initially this model did 
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not receive much attention, despite a successful demon- 
stration of this model for neutron scattering from various 
nuclei by Peterson 0, E] . There were some attempts to 
put this Ramsauer model on a sound theoretical basis 
lH IH 13 13 (see refer ences therein). The neutron to- 
tal cross sections have thus been well studied using this 
model, over a wide range of nuclear masses as well as neu- 
tron energies up to 500 MeV E E E E EE E . 
Various cross sections used in the Ramsauer model are 
given below. 



5000 



Otot 



2ir (R + X) 2 {l-acos(3) 



(5) 



= 7r(i? + A) 2 ((1 - a cos/?) 2 + {a sin/3) 2 (6) 



I. ANALYSIS OF TOTAL CROSS SECTIONS 

We performed the Ramsauer model fits to total cross 
sections by using Eq.(5,7) for a set of seven nuclei with 
a x 2 /N — 653.762 . The R, a, (3 are functions of atomic 
mass number and the center of mass energy. 

otot = 2vr (i? + A) 2 (l-acos/3) (7) 

a = «oM 1/3 ; 13/ AP = a ^b + c E - Ve) (8) 

r = 1.25519, a = 1.42042, p = 0.29361 (9) 
a = 3.0, b = 5.84633, c = 0.97676 

The a to t cross sections are shown in Fig. 1 with solid 
lines representing Ramsauer model fits using Eqs.(7-9) 
and the symbols represent the results from the optical 
model code SCAT2 |7J. The fits are obtained with total 
of six free parameters as given in Eq.(9), over wide mass 
range of 27 Al to 208 Pb and also covering the neutron 
energy region (E cm ) of 50-250 MeV. Similar Ramsauer 
model fits to total cross sections were already shown by 
various groups 0113 (and see the references therein). 
In these works, the normalised (relative to 208 Pb) cross 
sections were shown for various nuclei (I?) . Our func- 
tional dependence on energy and mass given in Eq.(8) 
with six global parameters was able to reproduce the 
SCAT2 results very well. The form of f3 given in Eq.(8) 
is chosen such that (3/A p is a global phase function for all 
nuclei, as shown versus energy in Figure 2. We have cho- 
sen this form of a, /3 which directly scale by some power 
of mass rather than a polynomial function of mass in or- 
der to obtain a global function. The energy dependence 
of global (3 function has also been well discussed in liter- 
ature. The parameter bo is called potential depth and cq 
(~ 0.90) is the non local parameter. The atomic number 
dependence of (3 can also be introduced through these 
parameters, however at the cost of global function. 
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FIG. 1: Ramsauer model fits (solid lines) to total scatter- 
ing cross sections of SCAT2 (symbols) versus E cm , using Eq. 
(7). The six parameters required are mentioned in the text 
in Eqs.(8,9). The curves from bottom to top are respectively 
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FIG. 2: The global f3 function required for Ramsauer model 
fits to total scattering cross sections using Eqs. (8,9) 



II. RELATIONS BETWEEN VARIOUS CROSS 
SECTIONS 

The total scattering cross section is defined as a sum of 
shape elastic and reaction cross sections given in Eq.(lO). 



CTtot = CTse + OVeac (10) 

Further, it has been observed [20I l2ll | earlier that the 
ratio (see Eq.(ll)) of zero degree elastic cross section 
and the shape elastic cross sections shows linearity with 
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FIG. 3: The 70 factor defined (Eq.(Il)) as the ratio of zero 
degree angular distributions and shape elastic distributions. 
This can be fitted by a functional form given in Eq.(12). 



energy. In the present work, we study this linearity in the 
energy range of 50-250 MeV. We observed that the linear 
relation is quite good at higher energy region above 150 
MeV and this linearity remains valid over a large energy 
range at high energies. These ratios from SCAT2 for 
various nuclei are shown in Fig. 3 by symbols. They 
have been parameterised well (see solid lines in figure) 
by a function of the form given by Eq.(12). Making use 
of the Ramsauer assumption for these cross sections in 
Eqs.(2,4), it can be shown that the ratio has a simple 
form as in Eq.(12). This is because, the ratio involves 
{imax + 1) where we use £ max = x + §• For To nts we 
used R = lAMA^fm in Eq.(12). 
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Therefore, inverting this equation and using the known 
shape elastic cross sections, we can obtain the zero degree 
elastic cross sections for a given energy. Therefore, 70 
serves as a relation between er se and <j{6 = 0). 
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III. ANALYSIS OF REACTION CROSS 
SECTIONS 

As mentioned in the introduction, we need to fit only 
two quantities and the other two will be derived. In sec- 
tion I, we parameterised the <Jtot(A, E) for all nuclei. The 
two important relations between various cross sections 
have been discussed in the previous section. Therefore, 
it suffices to parameterise one of the three remaining 
cross sections (7 reac ,a(8 = 0),a se . In this section we 
analyse the cr reoc , the reaction cross sections and obtain 
<t(9 = 0),a se as derived quantities. As shown in Eq.(3), 
the reaction cross sections involves only the transmission 
function i.e., (1-|%| 2 ). Unlike the cr tot ,cr se and the an- 
gular distributions, the reaction cross sections do not 
involve /3 phase function. Therefore, we proceed to obtain 
the empirical fits to the reaction cross sections predicted 
by SCAT2 code. For this fits, we used the functional form 
and parameters given below, yielding a x 2 /-^=323.641 



(Jreac = 7t (i£ + 1.5A) e~ a ; a = a VE/A p (14) 
a n = 0.20873 ; r = 1.38413 ; p = 0.42507 (15) 

It should be noted that the fitted parameters of Eq.(15) 
may not be equal to the corresponding parameters of 
Eq.(9), as these are effective parameters representing dif- 
ferent phenomenological functional forms. One should 
also note that the multiplying radius dependence of 
these two equations are respectively, 27r(i? + A) 2 and 
7r (R + 1.5A) . These cr reac fits using three parameters 
given in Eqs.(14,15), are shown in Fig. 4 for all the nu- 
clei studied. Figures (1,4) show that these cross sections 
of SCAT2 (symbols) are well reproduced by our empir- 
ical fitting functional forms and their parameters (lines 
in figures). Hence, we proceed to utilise these cr reac fits 
to obtain other cross sections predicted by SCAT2 code. 

The difference of the two parameterised cross sections 
gives shape elastic scattering as given by, a se =a to t — &reac 
as shown in Fig. 5(a). Making use of the 70 and the 
a se of Fig. 5(a), the angular distribution at zero degree 
versus energy are obtained as shown in Fig. 5(b). In 
summary, we have fitted total scattering cross sections 
and reaction excitation functions and hence we derived 
cr se , <Jeias {8 = 0) . These two derived cross sections also 
agree very well with the SCAT2 code predictions. 



IV. WICKS LIMIT MODIFICATION FOR 
FORWARD ELASTIC SCATTERING 

In this section, we describe the parameterization of en- 
ergy dependence of zero degree angular distributions and 
hence the a r eac,&se will be derived quantities. It is well 
known that the optical theorem relates the total scatter- 
ing cross section to the imaginary part of forward scatter- 
ing amplitude is given by Eq. (16). However, the forward 
elastic angular distributions are defined in Eq. (17). In 
the Wick's limit [22], the real part of the forward scat- 
tering amplitude is neglected as in Eq. (18). Therefore, 
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FIG. 4: Phenomenological fits (lines) to reaction cross sec- 
tions of optical model code SCAT2. The SCAT2 data (sym- 
bols) is fitted by exponential energy dependence multiplied 
by (R + 1.5X) 2 (see text). The fits need three parameters 
that are given in text. The curves from bottom to top are 
respectively for Al,Ca,Cu,Y,Sn,W,Pb. 
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the fractional deviation (5) of true angular distribution 
at zero degree for a given energy from the Wick's limit is 
defined in Eq. (19). 
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Whenever, 6 is small, Wick's limit gives a good approx- 
imation to zero degree angular distributions. In depth 
study of this deviation function has been performed for 
several nuclei by Dietrich et. at., |l9j . They have pre- 
sented the deviation function derived from optical model 
predictions and Wick's limit values. In the present work, 
following [19j , we have constructed the deviation function 
over 50-250 MeV center of mass energy range for seven 
target nuclei. We observed that these can be parame- 
terised in terms of a Woods-Saxon function of energy, 
whose parameters are nuclear mass number dependent. 
Given below are the functional form for the deviation 



FIG. 5: The derived cross sections for shape elastic (a) and 
the zero degree elastic angular distributions (b). The shape 
elastic is obtained as difference of total and reaction cross 
sections. The zero degree cross sections are obtained by using 
70 factor form shape elastic data. 



function and its parameters optimised for all nuclei. 
_ (V, + V 2 E) 

° ~ 1 + e (E -E)/a 

E Q = 61.557 + 0.47825A- 9.90931 1Q- 4 A 2 (21) 



V 2 = 0.00217- 2.26414 10 
Vi = 0.15 + 143.04/A 2 
a = 13.48384 
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A + 5.164 10~ 8 A 2 



The deviation function derived from SCAT2 and its em- 
pirical fits are shown in Fig. 6. As shown in figure, 
the function with nine free parameters listed above, gives 
very good fits to data of SCAT2. We have also tried a 
logarithmic function, which gives similar quality of fits to 
deviation function with fewer number of parameters. It 
should be noted here that the S function can be given by 
some empirical prescription for all nuclei, very much sim- 
ilar to the fits of total cross sections. Once, 5 of Eq.(19) 
is known, using a™ of Eq.(18), the zero degree elastic 
angular distribution versus energy can be obtained, as 
given by ctq = (1 + ^)< j q / ■ Consequent to knowing ctq, 
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FIG. 6: The empirical fits of deviation function versus energy, 
as explained in the text. 



we can use 70 of Eq.(ll) to get the shape elastic cross 
sections a se . By subtracting the o se from total cross 
section, we can fix the reaction cross section versus en- 
ergy. We followed this scheme and the results are shown 
in Fig. 7. The agreement between the derived quanti- 
ties from fits and the SCAT2 predictions is very good. 
It should be noted that the deviation function method 
accomplishes fitting of the zero degree angular distribu- 
tion without need for a phase (3 function. Otherwise, the 
angular distributions fits would need a phase function. 



V. ANALYSIS OF SHAPE ELASTIC CROSS 
SECTIONS 



In this section, we show the parameterization of the 
shape elastic cross sections, which need a phase (3 func- 
tion. The parameterization is performed very much sim- 
ilar to the fits of total cross sections. We followed Eq. 
(6) for the functional form and searched for the param- 
eters. The R, a, (3 are functions of atomic mass number 
and the center of mass energy. The best fit functional 
form and the parameters are given below which yielded 
X 2 /N — 1222.97 However, in these parameterizations, 
the p function is more complicated, while maintaining 
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FIG. 7: The derived cross sections for shape elastic and the 
reaction cross sections. 

the functional forms for other parameters. 
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= 7T (R+ 1.5A) 2 , r = 1.20815 
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= ao/A 1 / 3 ;a 2 =a 20 /A 2 / 3 
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P 


= A p a (V&o + coE - VI) 


(25) 




= 1.08627, a 20 = 12.70 6 32, p = 0.3146 


(26) 


fl 


= 0.5534, 6 = 33.0085, c = 0.8626 





The shape elastic cross sections fits are shown in Fig. 8 
for all the systems, and the quality of fits is good. Follow- 
ing the shape elastic parameterization, the forward an- 
gular distribution can be obtained using 70 factor. The 
reaction is given as difference of total and shape elastic 
cross sections. The two derived quantities, cr reac , oq agree 
well with SCAT2 data (not shown here). In summary, we 
considered three methods in sections III, IV, V where two 
of the four cross sections were parameterised and the re- 
maining two were derived, as shown in the table. 

Conclusion 



6 



3000 




i i i i i i i i i i i i i i i i i i i 

50 100 150 200 250 

E 

cm 



FIG. 8: Empirical fits of shape elastic cross sections. The 
fitting functions and parameters are given in Eqs. (22-26) 



Method 


Fitted quantities 


Derived quantities 


1. 


&tot, &reac 


0" se , G 


2. 


Otot , O 


G se , &reac 


3. 


""tot, ""se 


Greaa 



TABLE I: List of parameterised cross sections and the derived 
quantities in sections of the text. In above ao represents zero 
degree angular distributions. 



We have analysed the predictions of optical model code 
SCAT2 for neutron nucleus scattering using Koning De- 
laroche potentials. We perform the Ramsauer model pa- 
rameterization of total scattering cross sections and de- 
rive the systematics of the Ramsauer parameters. We 
have proposed an empirical fits to the reaction cross sec- 
tions of SCAT2 and using this, the shape elastic cross 
sections have been obtained as the difference of these 
two empirical results. Using the 70 factor which is the 
ratio of forward elastic scattering to the shape elastic 
scattering, we determine the forward elastic angular dis- 
tributions. Thus our empirical method is able to predict 
the SCAT2 results with a good accuracy for all these 
four quantities, namely er tot , cr reac , c se , a{9 — 0). Wc 
have studied the Wick's limit of the forward elastic cross 
sections from SCAT2. We parameterised the deviation 
function of these angular distributions from Wick's limit 
and using these as a alternative method, all the remain- 
ing cross sections have been obtained which compare well 
with SCAT2 results. The present method is about em- 
pirical estimates of the various cross sections, therefore 
it will be of practical importance wherever high energy 
neutron scattering cross sections are required for exam- 
ple for the accelerator driven sub critical assemblies. 
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